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When analysing a system consisting of both dark matter and dark energy, an often used practice 
in the literature is to neglect the perturbations in the dark energy component. However, it has 
recently been argued, through the use of numerical simulations, that one cannot do so. In this work 
we show that by neglecting such perturbations one is implicitly making a choice of gauge. As such, 
one no longer has the freedom to choose, for example, a gauge comoving with the dark matter - in 
fact doing so will give erroneous, gauge dependent results. We obtain results consistent with the 
numerical simulations by using the formalism of cosmological perturbation theory, and thus without 
resorting to involved numerical calculations. 
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I. INTRODUCTION 

Cosmological data sets have increased in both their 
number and quality in recent years and we now hold a 
wealth of data about our universe. This allows us to 
confront the theory with observations and thus obtain a 
knowledge of the intricacies of the universe and its evo- 
lution better than ever before. Current data allows us to 
conclude that the majority of the energy density of the 
universe today is in the so called 'dark sector'. That is, 
the universe is currently predominantly made up of dark 
energy - the mysterious force responsible for the accel- 
erated expansion of the universe at the current epoch - 
and dark matter [l|-|Ij]. 

The nature of dark energy is arguably one of the most 
important open questions in modern cosmology. An 
obvious 'candidate' for the dark energy component is the 
cosmological constant, A. In fact, current observations 
are consistent with ACDM, the standard cosmological 
model in which the dark energy is a cosmological 
constant [l[. However, in light of the 'cosmological 
constant problem,' that the observed value of A differs 
from that predicted from fundamental theories by over 
a hundred orders of magnitude H, the cosmological 
constant is perhaps a less favourable candidate for dark 
energy. Thus, much recent work has been focussed on 
a dynamical agent for the dark energy with a negative 
pressure which will cause the universe to accelerate in its 
expansion (see Ref. @ and references contained within). 

One of the most useful tools that theoretical cosmolo- 
gists have at their disposal is cosmological perturbation 
theory. The basic premise is simple: since the universe 
appears homogeneous and isotropic on large scales, one 
starts with a Friedmann-Robertson- Walker spacetime as 
a background solution, and then adds small inhomogc- 
neous perturbations on top. (For a selection of papers 



and reviews on the topic see, e.g., Ref. [71-11 lj). Since 
one perturbs the geometry of the background spacetime, 
through Einstein's field equations, this invokes perturba- 
tions on its matter content. For a universe whose dark 
energy component is pure cosmological constant only the 
dark matter will be perturbed, but in the case of dynam- 
ical dark energy there will be perturbations to both the 
dark matter and the dark energy components. However, 
an often used practice in the literature (e.g. Ref. [l2[ 
and see references in [l3j ) is to ignore the perturbations 
in the dark energy and simply use the well-known evolu- 
tion equation given by 



2HS C - 4irGp c 6 c = . 



(1.1) 
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where 8 C = 5p c /p c is the dark matter density contrast 
and an overdot denotes a time derivative. This issue was 
addressed in Ref. [HI], where Park et. al. showed, using 
numerical methods, that dark energy perturbations are 
not negligible and so, in general, cannot be ignored. 
In this paper we address a similar question. Using 
the formalism of cosmological perturbation theory, we 
show that ignoring the dark energy perturbations is 
in fact a gauge choice, and so if we choose to set the 
dark energy perturbations to zero then we no longer 
have the freedom to choose a gauge in which the dark 
matter evolution equation takes the form of Eq. (|1.1|) 
- we are still left with a system of coupled differential 
equations. Since this work focusses on the formalism 
of cosmological perturbation theory, and mathematics, 
and not involved numerical calculations, we make no 
numerical approximations in obtaining this result. 

The paper is structured as follows. In the next sec- 
tion we derive and present the governing equations for 
the dark energy and dark matter system under consider- 
ation. This is followed, in Section IIII1 by a brief recap 
of the gauge transformations of matter and metric per- 
turbations. Finally, in Section llVI we present our results 
and conclude. 
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II. GOVERNING EQUATIONS 

We consider here only scalar perturbations to a flat 
(K = 0) Friedmann-Robertson- Walker spacetime with 
the line element 

ds 2 = -(l + 2<j))dt 2 + 2a 2 B, i dtdx l 

+a 2 [(l - 2ip)5 i:j + 2E tij ]dx i dx j , (2.1) 

where a = a(t) is the scale factor, <p is the lapse function, 
ip is the dimensionless curvature perturbation and E and 
B make up the scalar shear, a, as 



cr = a 2 E - aB . 



(2.2) 



All perturbations are functions of space and time (e.g. 
<p = 4>(x l , t)). Throughout this paper, Latin indices i, j, k, 
take the value 1, 2 or 3, Greek indices /Lt, v, A, cover the 
full spacetime indices and we denote a partial derivative 
with a subscript comma. 

We assume that the dark matter and dark energy are 
non-interacting, and so energy momentum conservation 
for each fluid gives 1 



(2.4) 



Then, we obtain an evolution equation for each fluid from 
the energy (temporal) component of Eq. (|2.4D . 

5 Pa + 3H(6 Pa + 8P a ) 
V 2 

+ (p a + P a )-2(y a + v) = 3(p a + P a )ifi, (2.5) 

where the covariant velocity potential of each fluid is de- 
fined as 



V a = a(v a + B) , 



(2.6) 



v a is the scalar velocity potential of the ath fluid and 
H = a/a is the Hubble parameter. 

Considering the dark matter fluid and the dark energy 
scalar field, respectively, Eq. (|2.5|) then gives 



6 c = 3ij- —(a + V c ), 
a z 



(2.7) 



5<p + 3HSip + U 



.W a 2 



Sip 



(3^ 



-2U^, (2. 



where we have treated the field as a perfect fluid with 
energy density and pressure 



/v = \v 2 + u tV> , 



(2.9) 
(2.10) 



We have also used the background evolution equation for 
the dark energy scalar field 



(p + 3Hcp + U tV = 0, 



and the fact that 



y = -— 

<Po 



(2.11) 



(2.12) 



There is also a momentum conservation equation, coming 
from the spatial component of Eq. (|2.4[) . corresponding 
to each fluid 



V a - 3Hc 2 a V a + 



5P a 



P a 



0. 



(2.13) 



where c 2 = P a /p a . 

The Einstein field equations give [16[ 

V 2 

3H(ip + H6) T (ip + Ha)+4wG6p = 0, (2.14) 

a 2, 

i}> + H(/) + At:G{p + P)V = 0, (2.15) 

cr + Fcr-^ + V^O, (2.16) 

ij) + 3Htp + H<j)+ (3H 2 + 2H)(j) - A-kGSP = , (2.17) 

where the total matter quantities arc defined as the sum 
of the quantity for each fluid/field, i.e. 



dp = Sp v + Sp c , 
SP = SP V , 
{p + P)V = (p v + P v )V p + p B V 



(2.18) 
(2.19) 
(2.20) 



and we have used the fact that the dark matter is pres- 
surclcss, i.e. P c = = 5P C . Note that, by treating the 
scalar field as a fluid, we can write the energy density and 
pressure perturbation for the dark energy, respectively, as 

5p v = (p5(p - 4>tp 2 + U !V 5ip , (2.21) 
SP V = ipStp - <j)ip 2 - U iV S<p . (2.22) 

Introducing a new variable Z, both for notational con- 
venience and to assist with the following calculations, 
defined as 



1 This assumption is for brevity, and to allow us to deal with more 
manageable equations. The results highlighted in this work will 
still hold in the case of interacting fluids for which the overall 
energy-momentum tensor is covariantly conserved, but the com- 
ponents obey 

VrTf a)V = Q {a)v , (2.3) 
where <2( a )„ is the energy-momentum transfer to the ath fluid 



Z = 3(^ + H<f>) - —a , 



(2.23) 



we can rewrite Eq. ([2T7) , using Eqs. (|2TT4]) and ([2~T5]l . 

as 

/ • V 2 \ 

Z + 2HZ + (3ff + — J cj) = AttG(6p + 35P) . (2.24) 
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Then, from Eq. (f2~7) . 



such that the coordinates transform as 



Z = S c + 3H<f)+ —V c , 



(2.25) 



and from Eq. (|2.13p for the dark matter fluid (for which 



0), 



(2.26) 



Differentiating Eq. (|2.25|) gives 

Z = 5 C + 3(H<f>J+^-(v c - 2#K) ■ (2.27) 
Substituting this into Eq. (|2.24[) gives 

5 C + 2HS C - AttG Pc 5 c = 8irG(2(pSip - U lV 5tp) 

+ V C (6(H + H 2 ) + mTrGip 2 ) + 3HV C . (2.28) 



The evolution equation for the field is then obtained 
solely in terms of matter perturbations from Eq. (|2.8[) 
by using Eqs. (|2T2"6"|) and (pT2"5j) : 



Sip + 3HSp + ( U iVV ) Sip 



(2.29) 



V 



<p( S c + —V c -V c )+ 2U^V C 



It is worth restating that we have derived these equa- 
tions in a general form without fixing a gauge. If we set 
the dark matter velocity to zero, i.e. V c = 0, then they 
reduce to those presented in, for example, Refs. [HI, E3] 
(for the case of a zero energy- momentum transfer). 



III. GAUGE TRANSFORMATIONS 

As mentioned above, so far the equations have been 
presented without choosing a gauge. In order to choose a 
gauge, we need to first consider how the variables change 
under a gauge transformation, that is, a change of co- 
ordinates in the physical spacetime while holding the 
background coordinates fixed. Recall that, in general 
the gauge transformation of an arbitrary tensor field T 
is given by the exponential map [ill E3 



e^T. 



(3.1) 



This can then be expanded to give the gauge transfor- 
mation at linear order in perturbation theory as 



ST = ST + £^T , 



(3.2) 



where To is the value of the tensor field T in the back- 
ground and £^ denotes the Lie derivative with respect 
to the gauge transformation generating vector, This 
can be split up as 



t = t — a , 

x l = x l - /r 



(3.4) 
(3.5) 



Then, one can show that scalar quantities such as the 
field perturbation transform as 



Sip = Sip + (pa , 
the density contrast transforms as 

o c = <J C H a, 

Pc 



(3.6) 



(3.7) 



and the components of the velocity potential as 

v a = v a — a/3 . (3-8) 

Furthermore, by considering the transformation be- 
haviour of the metric tensor gives us the following trans- 
formation rules for the scalar metric perturbations 



(3-9) 
(3.10) 

(3.11) 

(3.12) 

(3.13) 



Finally Eq. p.8[) along with Eq. (|2.6[) . gives the trans- 
formation behaviour of the components of the covariant 
velocity potential 



ip = ip — Ha , 

B = B - -a + a/3 , 
a 

E = E + /3, 
and so the scalar shear transforms as 
a = a + a . 



V a = V a - a . 



(3.14) 



= (",/?/) 



(3-3) 



IV. DISCUSSION AND CONCLUSIONS 

When splitting the spacetime into a background and 
a perturbation, as is done in cosmological perturbation 
theory, one introduces a gauge problem. That is, there 
exist spurious gauge modes which must be removed in 
order to draw physically meaningful conclusions. This 
problem arises because, while general relativity is a fully 
covariant theory, this splitting is not a covariant process 
(see, e.g., Ref. [18||). In order to resolve this issue Bardeen 
proposed a systematic method of dealing with the gauge 
modes, by considering the transformation behaviour of 
quantities and constructing gauge invariant variables 
fl9j . We follow that approach here, by considering the 
gauge transformations of perturbations presented above 
in Section [TTT1 and making gauge choices such that we can 
remove the components of the generating vector a and 
(3. Any remaining quantities will then be gauge invariant. 
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Turning now to the case at hand, choosing the gauge 
in which the perturbation in the dark energy field is zero, 
Sip = 0, fixes a as 



8 ip 



(4.1) 



Since none of the gauge transformations of the quantities 
involved in the governing equations depend upon /?, we 
do not need to consider fixing this explicitly here. (Of 
course, one can rigorously fix /3 by choosing a suitable 
gauge condition, for example, by setting E = 0.) The 
governing equations in this gauge are then 

8 C + 2H8 C - 4irGp c S c 

= V C (&{H + H 2 ) + l&TrGf 2 ^ +3HV C (4.2) 

(4.3) 



V c . - —r^Vc - ^V c = S c , 



where the hat denotes that the variables are evaluated 
in the uniform field fluctuation gauge. That is, in this 
gauge, 8 C and V c are gauge invariant variables defined as 



S c = 5 C ^rSf , 

Pcf 



K = K + — . (4.4) 



Alternatively, choosing a gauge comoving with the 
dark matter, in which V c = fixes the generating vector 
as 

a = V c , (4.5) 
and reduces the governing equations to 



By studying the above systems of equations, it is evi- 
dent that choosing the dark energy field perturbation to 
be zero is a well defined choice of gauge, reducing the 
governing equations to Eqs. (|4.2[) and (|4.3|) . Then, hav- 
ing done so, we are no longer allowed the freedom to 
make another choice of gauge. Alternatively, choosing a 
gauge comoving with the dark matter uses up the gauge 
freedom, and so we are not permitted to neglect the per- 
turbation in the dark energy field. In fact, doing so will 
result in erroneous gauge dependent results. It is clear- 
est to see why this is the case by considering the set of 
governing equations. By making our choice of gauge we 
are left with a set of equations which is gauge invariant: 
that is, performing a gauge transformation will leave the 
set of equations unchanged. However, by neglecting the 
perturbation in the dark energy field after having chosen 
the gauge comoving with the dark matter amounts to 
setting the right hand side of Eq. (|4.6p to zero. This re- 
sulting equation will, in general, then no longer be gauge 
invariant. 

Thus, we conclude that the dark energy perturbation 
must be considered in a system containing a mixture of 
dark matter and dark energy. Our result is consistent 
with that of Ref. [Hj], though we have shown this by 
simply using the formalism of cosmological perturbation 
theory instead of relying on more involved numerical cal- 
culations. 



S c + 2H5 C - 4nGp c S c = 8nG(2if6(p - U yV Sf) , (4.6) 
V 2 \ - 

U.^ j)5tf = tp5 c , (4.7) 



5f + 3H8f + ( 



where the bar denotes variables in the comoving gauge 
and we have 



5f = Sf + ifV c 



(4.8) 
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